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We compute the strangeness-conserving double beta decay of S~ hyperons, which is the 
only hadronic system that can undergo such decays. We consider both, the lepton number 
conserving X~ — > £ + e~e~z/z/ {PfizzJ) an d the lepton number violating £~ — > S + e~e~ 
(PPno v ) modes. The branching ratios of these (3(3 decays are suppressed at the level of 
10 -30 considering a light neutrino scenario in the case of the (3(3^ channel. The dynamical 
origin of such low rates and their possible enhancements are briefly discussed. Given its 
simplicity those decays can be used also for the purposes of illustrating the main features 
of double beta decays. 
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1. Introduction 



Neutrinoless double beta (/3/3oi/) decays would occur if a mechanism allows the violation 
of the total lepton number L by two units. Their observation in experiments will provide 
unambiguous evidence for physics beyond the Standard Model (SM). At present, the first 
evidences for physics beyond the SM come from the flavor oscillations in the neutrino 
sector that required to explain the deficit observed in solar [ [l| and atmospheric neutrinos 
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[ ||. Flavor oscillations of neutrinos do not require a change of the total lepton number 
(namely |AL| = 0), allow us to conclude that neutrinos are massive, but do not establish 
whether they are Dirac or Majorana particles. If neutrinos are Majorana particles, f3f3o u 
decays have to occur at some level; their observation will establish the Majorana nature 
of neutrinos beyond any doubt. 

The strangeness-conserving (AS = 0) /3(3 0l/ decays have been unsuccessfully searched in 
nuclear transitions for several decades. An intense activity in the experimental and theo- 
retical fronts [ |3|, f|, ||, |6], 0, |8| witness the importance of such decays as a sensitive probe of 
physics beyond the SM. Some examples of extensions of the SM that can induce contribu- 
tions to PPov decays are right-handed weak couplings as well as the Higgs exchange [ [J, 
right-handed weak coupling involving heavy Majorana neutrinos [ |TtJ, massless Majoron 



emission [ f|, |TT], [12|, |H| [L4|, and R-parity breaking in the supersymmetric models [[if 



From a theoretical point of view, (3/3o u decays in nuclei are limited in precision due a wide 
range of model-dependent predictions for the nuclear wavefunctions. In the present work, 
we study the double beta (both the lepton number-conserving fifov an d lepton number- 
violating (3(3q v ) strangeness-conserving decays of the £~ hyperon. The £~ hyperon is a 
unique system in hadron physics that can undergo strangeness-conserving double beta 
decays as it will be explained below. The hadronic matrix elements necessary for such 
calculations are well known and, therefore, can exhibit the underlying mechanisms for 
double beta decays in a more clean way. 

As we have mentioned before, the determination of an upper bound for the effective 
neutrino mass in nuclear j3/3 0u decays is limited by model-dependent evaluations of the 
nuclear matrix elements. Some of the difficulties we encounter in those calculations are 
the following: 

1. The nucleus is a many body system with many degrees of freedom; in practice there 
is not a well defined rule to choose the most relevant components to describe an 
specific excitation; 

2. The Hilbert space where nuclear models are worked out have a huge dimension 
requiring a lot of time consuming computational work; and 



3. The multipole expansion for the f3(3 0u decay amplitude is rather complex making 
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theoretical expressions difficult to manipulate. 

Despite these limitations in theoretical inputs, the large sensitivity of present experi- 
ments have been able to set strong constrains on the so-called effective Majorana mass 
term, (m ee ) = J^U^m^, where m Ul denote neutrino mass eigenstates. By assuming that 
fifiov in nuclei are mediated by the exchange of light Majorana neutrinos, the experimental 



upper bound on (m ee ) is 0.2 eV [ [17 



Before we proceed with our calculation, it is interesting to take a look at other reactions 
that can provide information on the violation of lepton number. The upper limits available 
on these rare processes can be used to set upper limits on the matrix elements (m Q/ g) = 
J2U a iU/3im Ul , where a,j3 = e,fi,r. Thus, muon to positron conversion in nuclei + 
(A, Z) (A,Z - 2)e+ [ ! gives (m eAt ) < 17(82) GeV depending on the spin of the 



initial proton pair; the production of three muons in neutrino-nucleon scattering [ [L9 
leads to the upper limit (m^) < 10 4 GeV; this limit has been slightly improved at HERA 
through the reaction e + p — > vlfl^X [ |5D| , giving (m /iAt )^4x 10 3 GeV and also for first time 
limits on the (mi T ) (connected with the r-sector) were given; finally, the non-observation 
of heavy Majorana neutrinos at various colliders [ ^TJ can also be used to set limits on 
the effective Majorana mass. Similarly, some rare kaon decays are also useful to constrain 
lepton number violating interactions. For instance, present bounds on the branching ratio 
of the K + — > [i + fi + 7i~ decay [ ^] translates into the upper limit (m w ) < 4 x 10 4 MeV [ 



23]. Those bounds from collision experiments and rare kaon decays are several orders of 



magnitude above the limit < 4.4 eV, inferred by an analysis [ 24] that combines 



experimental constrains from atmospheric and solar neutrino oscillation and the tritium 



beta decay end-point experiment [ 25, 26 



Conversely, we can use the existing bounds from lepton number violating processes 
to set upper limits on the branching ratios of some rare kaon decay. Thus, the upper 
limit B(K + — > yU + /i + 7r~) ^ 10~ 30 (10~ 19 ) [ 2^] can be obtained in models with a light 
(heavy) neutrino scenario, while B(K + — > e + e + n~) ^ 1CT 26 [ 2S] can be derived from 



upper bounds from nuclear (5(3q v decays. These values are well below the sensitivities of 
current experiments; for example the limits reported by the E865 experiment [ ^3] are 
B(K + e + e + 7r-) < 6.4 x 1(T 10 , and B(K + /i+^+vr-) < 3 x 1(T 9 . This shows that 
we are far from detecting such processes, in spite of the special window in the hundred 
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MeV region [ , and that the sensitivity of such processes are well below the previously 
discussed (3(3q v case. Nevertheless it is important to pursue searches for | AL\ = 2 processes 
since they would lead to nonvanishing results even if nuclear (3(3 0v decay turns out to vanish 
or becomes extremely suppressed. 

In this Letter we study the two double beta decays of the S~ hyperon, namely: the 
lepton number conserving £~ — > £ + e~e _ z/z/ (/3/3e 2i/ ) an d the lepton number violating 
£~ — > £ + e~e~ (/3/3s 0l/ , AL = 2) channels. The Feynman graphs for the (3(3-^ 2v and 
PPz „ decays, indicating the two intermediate states, are shown in Figure la and lb, 
respectively. To our knowledge, these decay modes have not been reported previously in 
the literature and even an experimental upper limit is not available (an upper bound of 
order 10 -4 was reported in ref. [ for the branching ratio of the related strangeness- 
changing S~ — > pfi~[i~ double beta decay). The isotriplet of £ hyperons S°, £~) is 
a unique system of hadrons that can undergo strangeness-conserving double beta decays. 
Actually, the E~ and the E + are not antiparticles of each other; they have a mass splitting 



of m£- — = 8.08 MeV [ pD| , which allows a sufficient phase space for (3(3 decays. An 
interesting feature of this system is that we can identify only two well defined intermediate 
states (S° and A) that can give important contributions to the (3(3 transitions under 
consideration. One of these states (the S°) lies in the middle of the S + levels while 
the other (the A) is around 77 MeV below them. Namely, one of these intermediate states 
can be real while the other is virtual. 

On another hand, the hadronic matrix elements required for the calculations are domi- 
nated by a couple of axial and vector form factors, which have been computed by several 
authors in the literature [ [H], with good agreement among them. This is an important 
advantage over the evaluation of nuclear transition matrix elements. All these characteris- 
tics of the S hyperon system make very interesting the study of their double beta decays, 
even as a possible textbook example to learn the formalism of (3(3 transitions applied to 
free baryons. In the following we consider each case separately. 



2. Lepton number-conserving (3(3 decay 



One of the four Feynman graphs corresponding to the (3(3y^v decay is shown in Fig. la 
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(the other three contributions are obtained under proper antisymmetrization with respect 
to final state electrons and antineutrinos). 




E°, A 




v=v 



(a) 



(b) 



Figure 1. Lowest order diagrams contributing to the E decay for (a) two neutrino double 
beta decay mode, and (b) neutrinoless double beta decay mode. 



The effective four-fermion weak Hamiltonian acting at each vertex has the usual current- 
current form [ ||, |32 



G 

H w = -j^hf + h - c - 
where 



(1) 



(2) 



is the baryonic current operator underlying the B —+ B' transition and 



3n = i>el^ ~ 7s)Vv, 



(3) 



is the V — A leptonic current operatoiQ. In eq. (2), fi and g\ are dimensionless vector 
and axial-vector form factors for the B — > B' transition. In our approximation we have 
neglected their momentum transfer dependence and we have also neglected the small 
contributions of the induced magnetic and scalar form factors for the vector and axial 



1 We do not consider here the possible existence of right-handed neutrinos, since their contribution to the 
nuclear double beta decay is negligible [ H || . 
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currents. The effective weak coupling constant is G = GpV u d, where Gf = (1.16639 ± 
0.00003) x 10 -11 MeV~ 2 is the Fermi constant and V u d the relevant element of the Cabibbo- 
Kobayashi-Maskawa mixing matrix. 

The rate for the /3/3e 2 „ decay from the initial X~ = A to the final S + = B hyperon is 
given in the rest frame of A by (we use natural units, i.e., h = c = m e = 1) 

where uia (cb) is the energy of the initial (final) hyperon, and and and q«) denote 

the energy and momentum of the electron (neutrino). The decay amplitude reads 



M 2v = [1-P{e 1 e 2 ))[l-P{u 1 u 2 )) 

(Pb] e\e 2 v\v 2 \E w \r], Q] eit>i)(r7, Q; ejA\Hw\pA[ 
n i^ A m A ~ e Q (v) -ex-uJx 



with cq(?7) = yQ 2 + m v and Q = — pi — qi = Pb + P2 + <l2 being the energy and 
momentum of the intermediate baryon state, and P(x\x 2 ) the operator that exchanges 
xi with x 2 . Thus, after introducing eqs. (HD-(|D in eq. (|5|) we get 

M*, = -[1-PM)[1-P{v 1 » 2 )} V ^g^)Wff)^^gi.ft.^) |(6) 
2 m A - e Q {rj) - e x - u x 

where 

Bh(Pi,Pf\ x ) = u(Pf)iMxv + 9x v 75)u(pi) } (7) 
is the baryonic matrix element between states with momentum pi and pp, and 

L^(pi,qi,P2,q2) = Uefahni 1 -TsKfeK^ih^ 1 -75)wp(?i), (8) 

is the leptonic tensor. The values for the form factors = fi(Trrj), = ^(E - ^), 
/bjj = fiij]^) and = (7i(?7S + ) at zero momentum transfer are summarized in Table 

Given the small mass difference between the relevant hyperon states, we can use the 
non-relativistic impulse approximation for the baryonic current [ ^ . The small mass 
difference between the £ hyperon states is responsible of the suppression of the decay via 
the real intermediate S° particle, while the decay through an intermediate on-shell A state 
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Table 1 

Hyperon form factors at zero momentum transfer (see Table III on Ref. [ 32]). 



V 


fA V 


9a v 


/fir? 


9b v 


A 





-0.60 





-0.60 


S° 


1.41 


-0.69 


-1.41 


0.69 



is forbidden since m\ < mj. Keeping only the usually called Fermi and Gamow- Teller 
operators we have 

Bp{pi, Pf; X) = xln sp {fxr,9^ - gx V ^k9iJ,k)Xm ai , (9) 

where sj and sp denote the spin of the initial and final baryons in the I —> F transition. In 
the spirit of the non-relativistic approximation the energy denominator in eq. (§) can be 
also simplified making €Q(r]) = m v . Note that in the case of the E° intermediate state, eq. 
(6) exhibits a singularity when ei+cui = itla — m^o = 4.8 MeV (note that this singularity 
does not appear for the A intermediate state because uia — rn\ = 81.7 MeV e\ + lo\). 
This singularity can be cured by taking into account the finite width (T s o = 8.89 KeV [ 
|3~0|] ) of the E° intermediate state. Therefore, we will define: 

• In the case rj = S°: 

-> p— = hvpfa, Wi). (10) 

m A - m s o -e x -u)\ m A - m E o - e x - u x + i-^- 



In the case rj = A: 
1 

ITLa — 1TL\ — €\ — U)\ 



h K (ei,ui). (II) 



After a lengthy but straightforward calculation we can write the decay rate as follows: 

r ^ = SEE^i(W)iiW) • (12) 

07T' • -, , 
3=1 r\r]' 

The Cj(rjr]') terms denote quartic combinations of form factors, 

C\{t)T)') = f A-qf B-qf Ar)' f B-n' + ?>9 ArJ Br)9 Ar,' f Br)' + % f Ar}9 Br, f Ar,' 9 Br)> + 9 9 Ar,9 Br,9 Ar,' 9 Br,' ■, 
(^(w) = f Ar,f Br,f Ar,' f Br,' + ^9 Ar,f Br,f Ar,' 9Br,' + ^ f Ar,9 Br,9 Ar,' f Br,' ~ ^9 Ar,9Br,9 Ar,' 9Br,' ■, 
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^(W) — f Ar\j Bt)J At)' j Bt)> + 3 j Ujj j ? B V 9 'A v < 9 'B v > + 3 >j "at)' ] * Brf '9 'A v 9 'B v + "&9 ArJ Bt)J Ar)> 9Brf 

+ 3#Ar/ fsr! 1 fAr,gBr) + ?>9 Arif Br]9 Ar]< f Bf]< + 3J "~ Ari9Br)j ' Arf '9Br,' ~ ?>9 'Ar)9Br]9 'Ar? '9Br)> ', 

and 2j(r}r}') are the phase space factors defined as follows: 

/£0-l /-eo-ei reo-ei-€2 

Pidpi y i p 2 dp 2 J o 9i(eo - ei - e 2 - <?i) dgi 

x h v (e 2 ,e - ei - e 2 - qi)h*,(e 2 ,e - e x - e 2 - qi), 

^(VV) = / Pidpi P 2 dp2 gi(eo - ei - £2 - gi) rfgi 

Jo 



x i?e 



^(^2, e - ei - e 2 - <?i)/i*,(ei, gi) 



2s (W) = / p\dp\\ p\dp 2 \ ql{e Q - e x - e 2 - qxfdqi 

Ji Ji Jo 

h v (e 2 , e -e 1 -e 2 - q 1 )h*,(e 1 , e - e : - e 2 - g x ) 



re -«:i-e2 



x Re 



(13) 



(14) 



where e = rriA — m B . 

The numerical values for the factors entering in the expression of the decay rate, eq. 
(12), are given in Table 2. We can check that the main contribution comes from the term 
Ci(E°E )Ti(E E°), which includes the contribution of a real E° hyperon intermediate 
state. 



Table 2 

Numerical values for Cj(r]r]') and Xj(r}rj') (in units of MeV 9 ) for f3f3x 2u decay 



j 


C,-(AA) 


Jj(AA) 


^(EOE ) 


J^E ) 


^•(AE ) =C J 


;e°a) 


^•(AE ) +X j (E°A) 


1 


1.166 


5.46 x 10" 1 


11.672 


6.740 x 10 5 


-1.543 




4.106 x 10 


2 


-0.389 


5.46 x lO" 1 


8.952 


2.162 x 10 3 


0.514 




3.252 x 10 2 


3 


-0.389 


5.46 x 10" 1 


20.3101 


4.598 x 10 3 


-1.633 




4.170 x 10 



Using the values obtained in Table 2, we can compute the branching ratio from the rate 
in eq. (12). We obtain: 

B{(3fc 2 „) = 1.38 x 10~ 30 (1.36 x 10~ 30 ) . (15) 

Just for comparison we have shown within parenthesis the value corresponding to the 
contribution of the E° intermediate state. As expected, this contribution dominates al- 
most completely the decay rate. The branching ratio given above turns out to be very 
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suppressed due essentially to the large decay width of the S° hyperon appearing in eq. 
(10). As we know, the decay rate of the S° is 10 orders of magnitude larger than those 
of the charged £ hyperons because it can undergo the electromagnetic decay S° — > Aj. 



3. Neutrinoless double beta decay 

The decay rate for the (three-body) neutrinoless /3{3x 0v mode reads 
dT mu = it E \M 0u \ 2 5(m A - e B - ei - £2)7^77^7 • (16) 

The decay amplitude corresponding to the diagram in Figure lb (after proper anti- 
symmetrization with respect to identical electrons) is 



M 



[i-PM) EE/ 



d 4 q 

rj=T,°,A Su 



(Pb; eie 2 \H w \r}, Q(q)\exv){T}, Q(q);e 1 u\H w \p A ) 



(17) 



m A - ei - q - e Q (rj) 

The four-momentum of the intermediate state r\ is Q (q) = p A — Pi — q = Pb + P2 — q 



and its energy is eQ^rj) = ^/Q 2 (q) + m^. Introducing ([]]) into (0), and expanding the 
weak neutrino eigenstate as a mixture of light massive Majorana neutrino states, i.e., 
u„{q) = E« U e iu ui (q), we get 

M u = G 2 [l-P( ei e 2 )] (18) 

X E m ^^E^(^)^(/^ + 9Br,l5)I V (Plhv(fA V + 9A v l5)u(PA)L^ U (p 1 ,p 2 ), 

with the leptonic factor 

L»v(PuP2) = u e (pi)j^(l - ls)l v v%(p 2 ) ■ (19) 
The factor I v {pi) in eq. (|18D corresponds to the following loop integral (Q(q) = Pa—Pi—q)' 

I fa) = [Jil ^) + m v m 

Wl) J (2n)* (q* - mljmq) - m*y ^ 

which has a logarithmic divergence that we will manipulate in a simple cutoff procedure. 
After using Feynman parametrization techniques [ |33|, we get 

J,(Pi) = jl HWa - h)x + m v ) / o Ac {k2 \ dk M2r (21) 
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where A c is the cutoff energy and we have defined 



M 



m 2 (l — x) — (pa — Pi) 2 (l — x)x + m 2 V] x . 



(22) 



We stress here that the origin of this logarithmic divergence is related to the effective 
vertices we are using for the hadronic form factors. This divergence can in principle be 
cured by including the weak form factors which are expected to fall as 1/[(pa — Q(<?)] 2 
with q 2 — > oo in the dipole approximation, but their real behavior in the large q 2 limit are 
actually determined by QCD. In the case of the nuclear (3[3q v decays it is usual to assume 
that the divergence is driven by the internucleon distance, which sets A c ~ 1 GeV. This 
value can be identified with the largest momenta that the neutrino can carry which in 
turn is fixed by the lowest distance between two nucleons in the nucleus (A c ~ (2d) , 
where d is the nucleon radius). We will assume here that the average distance between 
the quarks within the hyperons are of the order of a typical hyperon (or nucleon) radius. 
Here, we adopt the point of view that our results are stable as far as they do not depend 
strongly on the specific value of the cutoff A c . 

The integration in eq. fl2T|) can be simplified by neglecting in eq. ( p2|) all lepton masses 
and momenta, which is consistent for the neutrino since we are assuming a light neutrino 
scenario. We obtain 



where 



(4tt) 5 



■fa)F v + m rj G rj 



(23) 



-(1 + m 2 )ln(m + m') — — (1 — m 2 )ln(l — m) — -m 2 ln(m) + — (1 — m 2 )ln(m') 



7T , 9N 2m' - (1 - mf (1 + m) 

i-(l + m 2 ) + — 

2 K J 2D 



1 — m\ fl+m 
arctg + arctg 



D 



D 



G v = -(1 + m)hi(m + m!) — (1 — m)ln(l — m) — mln(m) + -(1 — m)ln(m') 



Z7r(l + m) + 



2m' - (1 - m) 2 r 



D 



arctg 



1 — m 
D 



+ arctg 



1 + m 
D 



(24) 



In the above expressions, we have introduced the following dimensionless constants: m 



m 2 /m\, m! = K 2 /m 2 A and D = y4m' — (1 — m) 2 . 

In order to evaluate the unpolarized squared amplitude we will use the non-relativistic 
impulse approximation for the final baryon and take pb — 0) (we can not do the 
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same approximation for the intermediate baryon state in this case). The decay rate for 
the (3Po u transition becomes: 

TsO^ = ( m ee) -^f^Ov 
rm' 

(25) 

where (m ee ) = J^i^^U^ is the effective neutrino mass . The other factors appearing in 
eq. (|2"oD are defined as follows: 



^l(W) = (fAr)fBr) + gAngB^ifArffBr)' + gAn'gBr)'), 
^(VV) = (fArjfB^-gA^gB^ifAri'fBri'+gAq'gB-n'), 

^(vv) = {fAnfBr, - gA V gB V )(fA V 'fBn' - gAv'gBri'), (26) 
for the product of form factors, and 



Zo» = ] ^(^-^V^i- 1 )^- 6 !) 2 -!]^!. (27) 

for the phase space integral. 

The neutrinoless double beta decay rate depends on the cutoff A c and the neutrino 
effective mass (m ee ), which are free parameters in our model. Based on present bounds 
on electron neutrino mass ~ eV [ 17| , we show in Fig. the decay rate as a function 
of the cutoff A c , for a fixed neutrino mass of 10 eV. From this figure we can clearly see 
that our results are not very sensitive to the specific cutoff value in the region under 
consideration, which gives further support to the method employed for the logarithmic 
divergence. 

For illustrative purposes let us use A c = 1 GeV, which corresponds to a rough estimate 
of the inverse of the size of hyperons. The branching ratio in this case becomes 

B(f3f3 mu ) = 1.49 x 1(T 35 . (28) 

Thus, the neutrinoless double beta decay of the S~ hyperon is very suppressed and even 
smaller that the branching ratios of f3f3Q U decays of kaons. Note however that in the case 
of E hyperons we do not have a single neutrino (tree-level) intermediate state contribution 
as in the case of neutrinoless double beta kaon decays [ |27J . 
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Figure 2. Pflzois decay rate as a function of the cutoff A c , for (m ee ) = 10 eV. 

4. Summary and conclusions 

We have considered the strangeness-conserving double beta decays of £~ hyperons. 
Several characteristics make this particle a unique and interesting system to study (and 
to learn) such decays as an alternative to the corresponding decays in nuclei. First, the 
isotriplet of £ hyperons is the only system of hadrons that can undergo double beta decays 
since the J 3 = +1, —1 components have different masses. Second, we can identify a few 
intermediate states that dominate such decays where one state is in the middle and the 
other is below the initial and final hadronic levels. Finally, the uncertainties in hadronic 
matrix elements are much smaller than in the nuclear case. 

We study both the lepton number-conserving /3/3y;2v and the lepton number- violating 
decays. For the lepton number conserving decay we obtain B((3(3^ 2 u) = 1-38 x 1CT 30 
and for the neutrinoless decay we get B(PPxov) = 1-49 x 1CT 35 , standing in the light 
neutrino scenario in the last case. The suppression of the lepton number-conserving decay 
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is due to the large decay width of the S° intermediate state. Therefore, it would be nice 
to find an analogue system to the E hyperons (the baryons?) where the mechanism 
under consideration could produce an enhancement of the decay rate if the corresponding 
isotriplet is similar to the spectra of S hyperons. Our result for the neutrinoless double 
beta decay is almost insensitive to the specific value of the cutoff parameter used to 
regulate the divergent integrals. 

Our numerical results for the branching ratios may look discouraging. However, let us 
assume an hypothetical model where B{f3f3^ Qv ) = (10~ 20 ~ 10~ 25 )(m ee ) 2 ; in this very op- 
timistic scenario, an experimental upper limit of 10~ 8 would translate into the interesting 
bound (m ee ) < 1 ~ 300 MeV. On another hand, if there exists such a model that causes 
£~ to have faster neutrinoless double beta decays, the SM background due to the lepton 
number-conserving double beta decays would certainly be very small. 
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